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ABSTRACT 

The emergence of large-scale magnetic fields observed in the diffuse interstellar medium is 
explained by a turbulent dynamo. The underlying transport coefficients have previously been 
extracted from numerical simulations. So far, this was restricted to the kinematic regime, but 
we aim to extend our analysis into the realm of dynamically important fields. This marks an 
important step on which derived mean-field models rely to explain observed equipartition- 
strength fields. As in previous work, we diagnose turbulent transport coefficients by means 
of the test-field method. We derive quenching functions for the dynamo a effect, diamagnetic 
pumping, and turbulent diffusivity, which are compared with theoretical expectations. At late 
times, we observe the suppression of the vertical wind. Because this potentially affects the 
removal of small-scale magnetic helicity, new concerns arise about circumventing constraints 
imposed by the conservation of magnetic helicity at high magnetic Reynolds numbers. While 
present results cannot safely rule out this possibility, the issue only becomes important at late 
stages and is absent when the dynamo is quenched by the wind itself. 
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1 INTRODUCTION 

The presence of Galactic rotation and vertical stratification ren- 
ders interstellar turbul ence both inhomogeneo us and anisotropic. 
It is well known (e.g. iKrause & Radlerlll980l) that, in the frame- 
work of Reynolds-averaged equations, this will lead to the emer- 
gence of mean-field effects feeding on the spatial variations of 
the turbulen ce intensity and density field. Different th eoretical 
predictions dRiidiger & Kitchatinovlil993l : lFerrierej|l998h for the 
leading-order effects in the kine matic regime have pr eviously been 
tested by means of simulations dGressel et al .1120081) . We now ex- 
tend this approach into the domain of dynamically important mag- 
netic fields. The main aim of the present work is to extract the 
intricately non-linear dependence of the measured turbulence ef- 
fects on the strength of the Lorentz force - this is usually re- 
ferred to as "magnetic quenching". The knowledge of this func- 
tional dependence is required to evolve derived mean-field mod- 
els into saturation. This is important, because most observed spiral 
galaxies are found to harb our large-scale magnetic fields of near- 
equipartition strength (see I Fletcher! I2010] . for a topical compila- 
tion). Mean-field mod els are required to predict global symmetries 
of the magnetic field dBeck et alii 19961) . which is neither possible 
in local resolved simulations, nor in global simulations ignoring ef- 
fects occurring on u nresolved scales. There has been an active de- 
bate (see chapter 9 inlBrandenburg & Subramanian 2005[, hereafter 
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BS05) about severe constraints inflicted upo n the saturation level of 
the m ean-field dynamo, first addressed by I Vainshtein & Cattaneol 
dl992l) . These, potentially "catastrophic", constraints can arise as 
a result of the conservation of ma gnetic helicity a t large magnetic 
Reynolds numbers, R m. However, iKleeorin et al .1 d2002l) and later 
IShukurov et al. have demonstrated that such constraints can 

be alleviated by means of magnetic helicity fluxes. While such 
fluxes a re likely present in our simulations in the form of advec- 
tion (cf. lGressel et al.l2008l) . this may be restricted to the kinematic 
phase. For strong fields, we indeed observe the wind advection to be 
quenched itself in certain cases, and this will require careful atten- 
tion. In the present paper, we will however ignore this complication 
and focus on investigating the transition into the quenched regime. 
We will furthermore try to distinguish between different kinds of 
non-linearity that have been proposed. 



2 METHODS 



We employ the Nirvana-iii MHD code (Ziegler 2004) to perform 
direct simulations of the multi-phase interstellar medium (ISM) in 
a local patch (L v = L y = 0.8 kpc) of the Galaxy, but retaining the 
vertical structure of the disc (z e ±2.1 kpc). The grid resolution 
is A = 8.3 pc, an d we combine shearing-p eriodic boundary condi- 
tions in the plane dGressel & Zieglerl2007h with outflow at the ver- 
tical domain boundaries, where magnetic fields obey a normal-field 
condition. By means of localised injection of thermal energy, we 
mimic the effect of supernova explosions, thereby driving strongly 
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supersonic turbulence. Our numerical setup is similar to our previ- 
ous work but we here follow the evolution of the magnetic field into 
its saturation. During the course of the sim ulation, we emp l oy the 
test-field (TF) method first introduced by ISchrinner et alj d2005T) 
to infer turbulent closure parameters such as the a effect, turbu- 
lent diamagnetis m, and turbulent diffus ivity. The standard flavour 
of this method JSchrinner etafl l2007j)_ is now usually referred 
to as "quasi-kinematic" TF method iRheinhardt & Brandenburg! 
l20ld) . as it has been shown to remain fully valid into the non- 
kinema tic regime in the absenc e of magnetic background fluctu- 
ations l lBrandenburg et akll2008l) . Since our driving is dominated 
by kinetic forcing and we likely remain sub-cr itical to exciting a 
small-scale dynamo (see lMantere & Coldl2012h . this is warranted 
in our application. Whether, at higher Rm, the presence of a fluctu- 
ation dynamo will affect the evolution of the mean fields, remains 
open jBrandenburg. Sokoloff & Subramanianll2012h . 

We measure TF coefficients as a function of time, t, and verti- 
cal position, z. Due to the vertical stratification, the amplitude of the 
a effect significantly depends on position, as does the strength of 
the magnetic field. In other words, our system (unlike, e.g., forced 
turbulence in an unstratified triply-periodic box) has a rather poorly 
controlled background state and hence suffers from strong random 
fluctuations. Because of this, it is not straightforward to extract the 
quenching function. If we were to compile a naive scatter plot (as 
is usually done in a homogeneous setting) we are bound to confuse 
variations in a caused by fluctuations in the background gradients 
(i.e. in Vp and Vu 2 ras ) with the attenuation due to the non-linear 
quenching. To resolve this degeneracy, we have to make an assump- 
tion about the systematic variation with height. A very simple pos- 
sibility is to approximate the vertical profiles by a "power series" 
around z = 0, and for reasons of simplicity, we restrict ourselves to 
the lowest-order variation in z. For the a effect, and for the diamag- 
netic pumping, y z , we accordingly measure linear slopes within a 
defined central range in z. This is guided by inspection of Fig. [3] 
presented in the following section. The turbulent diffusivity, t]t, of 
course, has a dominant even parity with respect to z = 0, and we 
accordingly apply a quadratic fit. 

The vertical profiles of the magnetic and kinetic energy densi- 
ties are found to be roughly bell-shaped. As a measure for the rela- 
tive importance of the magnetic field, we introduce the quantity B = 
< 6/S cq ) as their ratio; here we have used B = |B| = (B 2 R + B 2 ) 1 ' 2 , 
and fi H] = (pop)' 12 Utms, with po the permeability of free space, and 
where overbars denote horizontal averages and angle brackets refer 
to rms values when averaging over (z, t). We remark that the weak 



vertical field, B-, remains constant throughout the run (a result of 
the periodic horizontal boundaries, conserving vertical flux), and 
hence does not influence the quenching function. 

Due to strong fluctuations within the data, and to obtain rea- 
sonable statistics, we bin the data into larger sections, for which 
we infer data points a(B), y z (B), and tjt(B) via linear regression 
as described above. Based on these data, we then attempt to 
fit the commonly adopted a lgebraic quenching formula (see e.g. 
Ilvanova & R uzmaikin 19771) 



with 3 independent parameters a , q„, and p„ in the case of a(B), 
and with analogous expressions for y z (B) and rjr{B), respectively. 
The following meanings can be ascribed to the three coefficients 
each: The subscript '0' refers to the unquenched amplitude, while 
the power p determines the abruptness of the transition into the 
quenched regime; q, taken to the power of —1/p, can be interpreted 
as the (approximate) relative amplitude B at which the quenching 
sets in. The latter observation points at a potential correlation be- 
tween these two parameters when fitting the data, and this is indeed 
confirmed for the case of a(B). 



3 RESULTS 

The model p r esente d in this work is identical to model "Q4" in 
iGressel et all J2008h . with the exception that it is initialised with 
a nano-Gauss vertical field rather than a horizontal oneQ The key 
parameters of the simulation are a supernova rate of a quarter the 
Galactic reference value of ctq = 30Myr~' kpc 2 , a rotation rate of 
four times 25 km s -1 kpc -1 , and background shear, q = — 1, equiva- 
lent to a flat rotation curve. F or further deta ils, we refer the reader 
to the extensive description in lGressell d2009l) . In the present model, 
the inclusion of a net-vertical flux was chosen to investigate a 
potential contribution from the magnetorot ational instability (see 
IPiontek & Ostrikedl2005l ; iKorpi et al.ll2010h . We here ignore this 
issue and focus on the quenching behaviour. The effect of different 
field topologies will be discussed in a forthcoming publication. 



1 Strictly speaking, B, is not a seed field but an imposed field. However, it 
induces rms fluctuations in Br, and B^ which then act as seed. 
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time t [Gyr] 

Figure 2. Time evolution of the fluctuating and regular magnetic field com- 
ponents, where overbars denote horizontal averages, and angle brackets re- 
fer to the rms variation within Br(z), and B^(z) at a given point in time. 
Note that the radial pitch angle (inset) is suppressed during saturation. 

3.1 Magnetic field evolution 

For our purposes, the only noticeable effect of the vertical seed field 
is that the growing mode takes a dipole-like symmetry with respect 
to the midplane, z = 0. This can be seen in Fig.Q] where we show 
the spatiotemporal structure of the horizontally averaged mean ra- 
dial field. Ultimately, the quadrupolar-like vertical symmetry will 
dominate, but not until the parity of the field has reversed twice. 
Figure Q] also reveals a remarkable feature caused by the multi- 
phase nature of the interstellar plasma: the occurrence of return 
flows, tracing-out the Galactic fountain. Due to shock compres- 
sion, the density is positively correlated with the magnetic field 
strength. As the dense clumps rain down into the Galactic poten- 
tial, the frozen-in field gets dragged along. A more comprehensive 
analysis of the correlation of mean magnetic fields with the ISM 
phases will require more sophisticated methods such as the ones 
recently proposed bv lGent et all l l2012h . 

If we apply a further rms average over the vertical direction, 
we obtain an estimate for the amplitude of the growing field, which 
is shown in Fig. [2] according to different constituents. The turbu- 
lent velocity tf lms reaches a statistically stationary state after less 
than lOOMyr, and with it the equipartition field ( S cq ). The fluctu- 
ating field ( B' ), and the radial and azimuthal mean fields show an 
extended phase of exponential growth, as indicated by the straight 
line. After the saturation of the mean-field dynamo (evident from 
( Bjt ) remaining constant), we observe that (B$) keeps growing due 
to stretching by the differential rotation. As is shown in the inset of 
Fig. [2] this has the consequence of a negligible radial pitch angle 
in the saturated state - which is at odds with observational evi- 
dence of s ignificant pi tch angles for equipartition-strength regular 
fields (see iBeckl 12013) . While this behaviour is characteristic for 
dynamos of the aSl type, it has been argued that the loss of pitch 
angle can be circumvente d if the dynamo is instead saturated by a 
wind l lElstner et ai]|2009f) . As we will substantiate shortly, this is 
not the case for the model "Q4" discussed here, but has in fact been 
observed at higher supernova rates (paper in preparation). 

3.2 Derivation of quenching functions 

We adopt the usual mean-field closure (cf . lBrandenburdr2005l) ap- 
plicable to the local box geometry, and assume that the turbulent 
electromotive force £(;, t) can be approximated by vertical profiles 
of the mean horizontal magnetic field B = (B R , S^) T , and its vertical 
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Figure 3. Time-averaged vertical profiles of the governing mean-field co- 
efficients in the kinematic regime (light colours, t e [0.1,0.3] Gyr), and in 
the fully quenched phase (dark colours, ( e [2.3, 2.5] Gyr). Panels show, 
(a) the helical dynamo effect, a^(z), (b) the diamagnetic turbulent pump- 
ing, y z (z) = 0.5 (a^R-aR^), and (c) the turbulent diffusivity, r]M(z). We also 
overlay the mean vertical flow (dotted), mirrored for the sake of easier com- 
parison. Note that, in both regimes, y- closely matches - U z in the dynamo- 
active region, rendering saturation via the wind mechanism unlikely. This 
is in-line with the loss of pitch angle seen in Fig. [2] 

gradients, i.e., 

£, = a-,j Bj - j]jjEjkid k Bi, i,je{R,(p},k = z- (2) 

We here focus on the relevant components of the tensors try and 
rjjj, namely the dominant contribution cr^(z) to the aQ dynamo, the 
vertical diamagnetic effect y z (z) = (a^g - a R ^,)/2, and the turbulent 
diffusivity t]^. These quantities are plotted in the three panels of 
Fig. [3] respectively. 

For each of the three effects, we show profiles at an early stage 
(light colours), representing the (magnetically unquenched) kine- 
matic regime, and at a late time, representing the fully quenched 
regime. We observe that all three profiles are suppressed in the 
central region, corresponding to the belt of strong fields seen in 
Fig.Q] We furthermore find the mean flow U : to be shut-off by the 
emerging magnetic field, explaining the absence of quenching via 
the wind, and hence the low pitch angle. This is somewhat unex- 
pected and will require careful further inspection - especially since 
the situation appears to change at higher supernova rate. 

The central aim of the current paper is to extract quenching 
functions. We here restrict ourselves to run "Q4" as it presents 
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Figure 4. Amplitude of the a effect as a function of relative field strength 
y3 = (B/B cq ). Note that quenching occurs significantly below equipartition. 
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Figure 5. Same as Fig. [4] but for the amplitude |y-| of the diamagnetic turbu- 
lent pumping term. Compared to oa^, the pumping appears to be quenched 
at somewhat higher relative field strength. 



us with the widest range of magnetic field strengths. Comparison 
with other runs, however, indicates that the derived relations should 
be fairly universal. To disentangle the vertical variation from the 
magnetic field dependence, we make a model assumption on the 
shape of the profiles (see Sect. [2] for details). For a^(z) we mea- 
sure the profile slope within |z| < 1.33 kpc; the resulting ampli- 
tude is plotted in Fig. [4] as a function of the relative field strength. 
Strictly speaking, the ordinate should have additional units kpc -1 , 
but we drop these for the sake of simplicity (and implicitly inter- 
pret the values as integral value at 1 kpc). With p a = 1.98 + 0.3, 
the fitted algebraic quenching function, cf. equation Q}, is found 
in close agreement with the co mmonly assumed quadratic (e.g. 
iRogachevskii & Kleeorir] |2000h depen dence on ft; bu t note that 
the classic result of Moffattl < 19721) and lRiidigerl J 19741) suggests a 
steeper scaling with p„ = 3 - notably for fields exceed ing equipar- 
tition, however. In the limit of low Rm. lSur et al. have found 
p = 2 for stochastic forcing, and p = 3 for steady forcing, which 
might provide an explanation. We, moreover, find a rather large co- 
efficient q a = 27 ± 14. This is in fact cause for some concern as 
one would naively expect q a ~ 10, based on the scale separation 
Lg/h — 1 kpc/0.1 kpc, as can be derived in the framework of dy- 
namical quenching (see sect. 9.3.1 in BS05). A value significantly 
larger than this would be indicative of "catastrophic" quenching, for 
which one can estimate q a ~ Rm. In previous work, we adopted a 
definition of Rm based on the global box scale, arriving at a (seem- 
ingly large) value of Rm' = L 2 Cl/q = 10, 000. In the current con- 
text, it is however more realistic to base Rm on the outer scale 
of the turbulence, i.e. k f ai 2n/l , and a typical rms velocity, e.g. 
Ums - 25-50 km s _1 . With this new definition, we obtain a value of 
Rm = u tms (k { qY l =! 75-125 which is still comfortably larger than 
the factor appearing in the quenching expression. Another com- 
monly used definition of Rm in this context is Rmr = qj/q. With 
r]i — 0.1 - 2.0 kpc km s~' in the region of interest, this makes for 
a poor indicator as we obtain Rm T a 15-300, with uncertain con- 
clusions. Nevertheless, the presence of catastrophic quenching can- 
not be safely ruled out here. Following the notion that alleviating 
helicity const raints requires a flux of helicity out of the dynamo- 
active region jKleeorin et alj2002h , its occurrence would, however, 
be consistent with the suppression of the wind profile (see dotted 
line in Fig.O at late times. Since the wind may be the deciding 
factor it becomes mandatory to better understand its exact physical 
origin. 
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Figure 6. Same as Fig. [4] but for the turbulent diffusivity 



For the diamagnetic term y z (z), which is more centrally con- 
fined, we infer slopes within \z\ < 0.5 kpc. Results are plotted in 
Fig. |U and with q y = 9.8 ± 1.6, and p y = 2.16 ± 0.25, they 
again agree with the assumption of a quadratic dependence ~ 1 . 
Like for or, we are at odds with the theoretical expectations of 
iKitchatinov &Rudigen < fl992h . who predict the same /? 3 scaling for 
y z - but recall the result by ISur et al. I <2007b . mentioned above, 
which might resolve this. Because diamagnetic pumping is footed 
on inhomogeneity of the turbulence, but not on lack of mirror sym- 
metry, it is not affected by issues of helicity conservation. This is 
consistent with a low value of q y . We moreover performed a cor- 
responding fit for the mean flow V-ifi) and obtain comparable pa- 
rameters. This is expected from the close correspondence between 
the two quantities seen in Fig. [3] 

A potential buoyant contribution to the pumping 
dKitchatinov & Riidigej 1 19921) appears unlikely in our case 
as the Galaxy (unlike the solar convection zone) is only moderately 
stratified in density. Similarly, tu rbulent transport a rising from 
small-scale magnetic fluctuations l lRadler et al]|2003h may only 
become important at higher Rm, due to small-scale dynamo ac- 
tion. Because we find no direct evidence for the latter, we attribute 
the observed attenuation to the quenching of the kinematic effect 
y z ~ -l/2Wq T . 

To conclude the results section, we present, in Fig. [6] the 
quenching function for the turbulent diffusivity qj (fl) . We note that 
due to background shear, the q tensor becomes anisotropic. We 



2 As a potential caveat, we want to point out that (as speculated upon ear- 
lier) the p and q parameters are found to be correlated at the 90% level. 



3 Whether this is expected to counteract the kinematic effect is, however, 
dependent on rotation (cf. sect. 10.3 in BS05). 
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however cannot address this issue with the current data. For practi- 
cal purposes, we restrict our analysis to rj^,- Unlike for the two pre- 
ceding quantities related to the a tensor, we now find a scaling ex- 
ponent closer to one, i.e. p 1} = 0.96±0.06. This is in good agree ment 
with p n = 1, as predicted bv lRogachevskii & Kleeorin! feOOd). and 
coinci des with the scaling suggestecHbv lBlackman & Brandenburg! 

Our result also agrees with the effective diffusivity (in 
the one-dimensional ca s e) quoted in t he lin e after eqn. (46) of 
iKitchatinov et alj Jl994h . lYousef et all ( l2003h have derived diffu- 
sivity quenching from simulations of forced turbulence. Their re- 
sults appear consistent with an assumed p,, = 2, based on which 
they state a best-fit value of q n = 8. 

4 DISCUSSION & CONCLUSIONS 

We have provided magnetic quenching functions for the attenuation 
of turbulent transport coefficients in supersonic interstellar turbu- 
lence. Our results are based on realistic simulations of the stratified 
ISM - with the excepti ons of neglecting a potent ial contribution 
from cosmic rays (CRs. lHanasz et "aH l2004l2006h . Nevertheless, 
the obtained quenching laws represent the first quantitative result 
derived from such a complex scenario and will be highly useful in 
the global modelling of the Galactic dynamo. 

Generally, the obtained results appear consistent with the pic- 
ture of purely kinematic effects, i.e., those arising out of correla- 
tions in the fluctuations of M rms . These correlations are then atten- 
uated by the presence of a strong mean field. A key question to 
address in future studies is, in how far correlations in background 
magnetic fluctuations (stemming, e.g., from a small-scale dynamo) 
will be of importance. Such contributions have usually been ig- 
nored since they are hard to capture. On the other hand, recent stud- 
ies have found that the small-scale dynamo is very hard to excite in 
the multi-phase ISM, which in turn challenges its importance. 

Independent of these, a potential influence of additional con- 
straints (arising from magnetic helicity conservation at high mag- 
netic Reynolds numbers) cannot safely be excluded at the current 
stage. Given the considerable uncertainty in determining the related 
q parameter, we conclude that settling this issue will require more 
detailed investigation of the scaling of our result with Rm, demand- 
ing for potentially very expensive co mputations. It is interes ting to 
note that, in comparable simulations. iDel Sordo et al" found 
the influence of a wind enhanced at higher Rm, pointing to a possi- 
ble solution. As far as our own simulations are concerned, we will 
need a thorough understanding of the conditions under which the 
vertical wind is attenuated as it may play a crucial role in alleviating 
helicity constraints. Because the inclusion of a CR component may 
potentially resolve this deficiency, we are nevertheless confident 
that the paradigm of a turbulent mean-field dynamo will prevail in 
explaining large-scale Galactic magnetic fields. 
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